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Nonequilibrium stationary states of thermodynamic systems dissipate a positive amount of energy 
per unit of time. If we consider transformations of such states that are realized by letting the driving 
depend on time, the amount of energy dissipated in an unbounded time window becomes then 
infinite. Following the general proposal by Oono and Paniconi and using results of the macroscopic 
fluctuation theory, we give a natural definition of a renormalized work performed along any given 
transformation. We then show that the renormalized work satisfies a Clausius inequality and prove 
that equality is achieved for very slow transformations, that is in the quasi static limit. We finally 
connect the renormalized work to the quasi potential of the macroscopic fluctuation theory that 
gives the probability of fluctuations in the stationary nonequilibrium ensemble. 
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A main goal of nonequilibrium thermodynamics is 
to construct analogues of thermodynamic potentials 
for nonequilibrium stationary states. These potentials 
should describe the typical macroscopic behavior of the 
system as well as the asymptotic probability of fluctu- 
ations. As it has been shown in this program can 
be implemented without the explicit knowledge of the 
stationary ensemble and requires as input the macro- 
scopic dynamical behavior of systems which can be char- 
acterized by the transport coefficients. This theory, now 
known as macroscopic fluctuation theory, is based on an 
extension of Einstein equilibrium fluctuation theory to 
stationary nonequilibrium states combined with a dy- 
namical point of view. It has been very powerful in 
studying concrete microscopic models but can be used 
also as a phenomenological theory. It has led to several 
new interesting predictions 

In this letter we develop a theory of thermodynamic 
transformations for nonequilibrium stationary states in 
the concrete case of driven diffusive systems. We thus 
consider an open system in contact with boundary reser- 
voirs and/or under the action of an external field and 
we analyze the situation in which the reservoirs and field 
vary with time driving the system from a nonequilibrium 
state to another one. For simplicity of notation, we re- 
strict to the case of a single conservation law, e.g., the 
conservation of the mass. The system is in contact with 
boundary reservoirs, characterized by their chemical po- 
tential A, and under the action of an external field E. 
We denote by A C M. d the bounded region occupied by 
the system, by x the macroscopic space coordinates and 
by t the macroscopic time. With respect to our previous 
work [J, [H , we consider the case in which A and E 
can depend explicitly on the time t. The macroscopic 
dynamics is given by the hydrodynamic equation for the 



density which satisfies the following general assumption, 
based on the notion of local equilibrium. For stochastic 
lattice gases these assumptions can be proven rigorously 
and the macroscopic transport coefficient can be charac- 
terized in terms of the underlying microscopic dynamics 

The macroscopic state is completely described by the 
local density u(t, x) and the associated current j(t, x). In 
the sequel we drop the dependence on the space coordi- 
nate x from the notation. The macroscopic evolution is 
given by the continuity equation 

d t u(t) + V-j(t) = 0, (1) 

together with the constitutive equation j{t) = J(t,u(t)) 
expressing the local current in function of the local den- 
sity. For driven diffusive systems the constitutive equa- 
tion takes the form 

J(t,p) = -D(p)Vp + X (p)E(t) (2) 

where the diffusion coefficient D(p) and the mobility x(p) 
are d x d positive matrices. In the case of time indepen- 
dent driving the right hand side does not depend explic- 
itly on time and we denote the current simply by J{p). 
The transport coefficients in ([2|) satisfy the local Einstein 
relation D{p) — /"(/?), where / is the equilibrium 
free energy per unit of volume. The interaction with the 
external reservoirs specify the boundary conditions for 
the evolution defined by {J)-©. Recalling that \{t) is 
the chemical potential of the reservoirs, this boundary 
condition reads f'(u(t,x)) = X(t,x), x e <9A. We shall 
also assume that when A and E do not depend on time 
there is a unique and globally attractive stationary so- 
lution for the flow defined by (P)-© that is denoted by 
P — P\,E- In particular, p~\ t E is the typical density profile 
in the stationary nonequilibrium state corresponding to 
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time independent chemical potential A and external field 
E. 

Fix time dependent paths X(t) of the chemical poten- 
tial and E(t) of the driving field. Given a density profile 
p, denote by (u(t),j(t)), t > 0, the solution of ©-© 
with initial condition p. Let WVt] = W[o,r] E, p) be 
the energy exchanged between the system and the exter- 
nal driving in the time interval [0,T], that is 



TT r 



[0,T] 



dxj(t)-E(t)}, 



da(x) A(t) j(t) ■ h 

(3) 



where h is the outer normal to A and da is the surface 
measure on OA. The first term on the right hand side is 
the energy provided by the reservoirs while the second is 
the energy provided by the external field. 

In view of the boundary conditions and the Einstein 
relation, by using the divergence theorem in ([3]), we de- 
duce that 

W [0tT] = F(u(T)) - F(p) 

f T f 1 ( 4 ) 

+ dt dxj{t)-x{u{t))- x j{t), 

JO J A 

where F is the equilibrium free energy functional, 

F(p)= [ dxf(p(x)). (5) 

J A 

Consider two stationary states corresponding to (time- 
independent) (Ao,-^) and (\i,E\) and denote by po — 
P\ ,E an d Pi — PXi.Et the associated density profiles. 
Such states can be either equilibrium or nonequilibrium 
states. We can drive the system from the initial state 
po at time t = to the final state p\ at time t = +oo 
by considering a time dependent forcing (X(t),E(t)) sat- 
isfying (A(0),£(0)) = {Xo,E ) and (A(+oo),£(+oo)) = 
(Ai,-Ei). As the second term on the right hand side of (|4j) 
is positive, by letting W = W^o.+oo) be the total energy 
exchanged in the transformation, we deduce the Clausius 
inequality 



W > AF — F(p~i) - F(po). 



(6) 



When the initial and final states are equilibrium states, 
e.g., the external field E vanishes and the chemical poten- 
tial A is constant, the inequality ((6]) is a standard formu- 
lation of the second law of thermodynamics. Moreover, 
by considering a sequence of transformations in which 
the variation of the driving becomes very slow, it is not 
difficult to show that equality in ([6]) is achieved in the 
quasi static limit, we refer to [10| for the details. On the 
other hand, in the present situation of nonequilibrium 
states the inequality (|6]) does not carry any information. 
Indeed, as nonequilibrium states support a non vanish- 
ing current, in the limit T — > +oo the second term on 



the right hand side of ((4]) becomes infinite so that the 
left hand side of (j6)) is infinite while the right hand side 
is bounded. By interpreting the ideas in [ll|, further de- 
veloped in 12, 13| . we next define a renormalized work 
for a which a significant Clausius inequality can be ob- 
tained also for nonequilibrium stationary states. 

To this aim, we recall the quasi potential, which is the 
key notion of the macroscopic fluctuation theory. Con- 
sider a system with time independent driving and let 
(u(t),j(t)) 7 t e [Ti,T 2 ] be a pair density-current satis- 
fying the continuity equation dfU + V ■ ] = 0. According 
to the basic principles of the macroscopic fluctuation the- 
ory 0j [3 j the probability of observing this path is given, 
up to a prefactor, by 



exp • 



£ d 0I[T u T, 



(7) 



where e is the scaling parameter, i.e., the ratio between 
the microscopic length scale (say the typical intermolec- 
ular distance) and the macroscopic one, /3 = 1/ kT (here 
T is the temperature), and the action functional I has 
the form 

imM^j) = \J t 2<it J A dx [x*) - J (m)] (8) 
■xim-'im-Jim}- 

In particular, if (u, j) solves ([I])-© then I[T lt T 2 ] 3) = 0- 
The above statement therefore implies that the typical 
behavior of the system is described by the hydrodynamic 
equations. The quasi potential is the functional on the 
set of density profiles defined by the variational problem 



V(p) = inf U{-oofi]{^,3) , w(0) = pj 



(9) 



where the infimum is carried out over all the trajectories 
satisfying the prescribed boundary condition. Namely, 
V(p) is the minimal action to bring the system from the 
typical density profile p to the fluctuation p. The proba- 
bility of a density profile p in the stationary nonequilib- 
rium ensemble is then given, up to a prefactor, by 



exp 



{-e- d pV{p)}. 



(10) 



In particular, the minimizer of V is the typical density 
profile p. For equilibrium states it can be shown [j| that 
V coincides, apart an affine transformation, with the free 
energy functional ([5]). Moreover, as shown in [l|, the 
functional V solves the stationary Hamilton- Jacobi equa- 
tion 

SV 6V f SV 

A op dp J A dp 

where SV/5p vanishes at the boundary dA and p satisfies 
the boundary condition f'(p(x)) = X(x), x S dA. 
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In the case of time independent driving, the current 
J(p) in ([2]) can be decomposed as J(p) = Js(p) + Ja(p), 
where 



Mp) = -x(p)v 



sv 

bp 



(12) 



and Ja(p) — J(p) ~ Js(p)- In view of the stationary 
Hamilton- Jacobi equation (|11|) , the above decomposition 
is orthogonal in the sense that for each density profile p 



dx J s {p) ■ x(p) 1 Ja(p) = 0. 



(13) 



We shall refer to Js(p) as the symmetric current and 
to J A {p) as the antisymmetric current. This terminol- 
ogy refers to symmetric and antisymmetric part of the 
microscopic dynamics [l|, 0, [l4| ■ We remark that J$ is 
proportional to the thermodynamic force and the above 
decomposition depends on the external driving. 

Since the quasi potential V is minimal in the stationary 
profile, we deduce that Js(p) — 0; namely, the stationary 
current is purely antisymmetric. In particular, J A {p) is 
the typical current in the stationary nonequilibrium en- 
semble associated and it is therefore experimentally ac- 
cessible. In view of the general formula (j4j for the total 
work, the amount of energy per unit of time needed to 
maintain the system in the stationary profile p is 



dx J A (p) ■ x(p) 1 Ja(p)- 



(14) 



We shall next consider time dependent driving and 
define a renormalized work by subtracting to the total 
work exchanged between the driving and the system the 
energy needed to maintain the stationary state. Fix, 
therefore, T > 0, a density profile p, and space-time de- 
pendent chemical potentials X(t) and external field E(t), 
t G [0,T], and let (u(t),j(t)), be the corresponding solu- 
tion of ([I])-© with initial condition p. Recalling ([13)) . 

WXU\,E,p) 



we define the renormalized work Wj™^ 



performed by the reservoirs and the external field in the 
time interval [0, T] as 

W( cn T] = W [0tT] 

rT f (15) 

dt / dxJ A (t,u(t)) ■ x(u(t))~ l J A (t, u(t)). 
la J A 

In this formula WVt] = W[o,T] (A, E, p) is given in ([3]), 

J(t,p) = J s {t,p) + J A (t,p), 

,8V X ( t ),E(t) (p) 



Mt,p) = -x(p)v- 



5p 



in which J(t,p) is given by ((2|) and V\u),E(t) is the quasi 
potential relative to the state X(t),E(t) with frozen t. 
Observe that the definition of the renormalized work in- 
volves the antisymmetric current J A (t) computed not at 



density profile p\(t),E{t) but at the solution u(t) of the 
time dependent hydrodynamic equation. That is, at time 
t we subtract the power the system would dissipate if its 
actual state u(t) were a stationary state. This choice, 
which is certainly reasonable for slow transformations, 
leads to a Clausius inequality. Indeed, by using ^ and 
the orthogonality between the symmetric and the anti- 
symmetric part of the current, 

W{ ™ T] (X,E,p) = F(u(T))-F(p) 

rT 

dt I dx J S (*,«(*)) • -*3 (*, "(*))■ 



Consider a density profile p and a space-time depen- 
dent chemical potential and external field (X(t), E(t)), 
t > 0, converging to (Ai,i?i) as t — > +oo. Let p\ = 
PXi.Et be the stationary profile associated to (Ai, E\) and 
{u{t),j(t)), t>0,be the solution of ©h© with initial 
condition p. Since u(T) converges to p\, the symmetric 
part of the current, Js(u(T)), relaxes as T — ► +oo to 
Js(pi) = 0. Under suitable assumptions on the trans- 
formation, the last integral in the previous formula is 
convergent as T — > +oo and we get 

W rcn (X, E, p) = F(pi) - F(p) 



dt / dx J s (*,«(t))-x(«(*)) _1 Js(t,tt(t)) 



(16) 



where F is the equilibrium free energy functional ([S]) . In 
particular, 



W mn (X,E,p) > F(pi) — F(p), 



(17) 



which is a meaningful version of the Clausius inequality 
for nonequilibrium states. Indeed, by considering a se- 
quence of transformations (X(t), E(t)) which vary on a 
time scale 1/5, we realize that the integrand on second 
term in the right hand side of (|T6|) is of order 5 2 while 
the integral essentially extends, due to the finite relax- 
ation time of the system, over an interval of order <5 _1 . 
Therefore in quasi static limit 5—^0 equality in (|17p is 
achieved. We refer to [lOj for more details. 

For special transformations, the integral in (|16p . which 
represents the excess work over a quasi static transfor- 
mations, can be related to the quasi potential. Consider 
at time t — a stationary nonequilibrium profile po cor- 
responding to some driving (Xo,Eq). The system is put 
in contact with new reservoirs at chemical potential Ai 
and a new external field E\ . For t > the system evolves 
according to the hydrodynamic equation ([I])-© with ini- 
tial condition p~o, time independent boundary condition 
Ai and external field E\. In particular, as t — > +oo the 
system relaxes to p\, the stationary density profile cor- 
responding to (X\,E\). A simple calculation shows that 
along such a path 



dt \ dx J s (t,u{t)) ■ xH^r 1 Js(t,u{t)) 
a Ja (lo) 

= V\ uEl (po) -V\ 1 ,e 1 (pi) = V\ uEl (po)- 
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This relationship shows that the quasi potential repre- 
sents the extra work necessary to reach the initial state 
from the final stationary state or the maximal useful work 
obtainable from the relaxation process. In equilibrium 
thermodynamics this quantity is known as availability 
(lif . Therefore, while a definition of thermodynamic po- 
tentials, that is functionals of the state of the system, 
does not appear possible in nonequilibrium thermody- 
namics, the quasi potential is the natural extension of 
the availability which is instead a function of the state of 
the system and of the environment. We refer to [loj for 
an interesting connection of the quasi potential with the 
relative entropy between nonequilibrium states. 

The definition of renormalized work we have intro- 
duced is natural within our approach and insures, as 
we have discussed, both its finiteness and the validity 
of Clausius inequality. It is however not simple to obtain 
as it requires the preliminary calculation of the quasi- 
potential, a generically nonlocal quantity. One may ask 
whether there exist other simpler ways of renormalizing 
in particular by subtracting a local counter term. To be 
concrete let us consider again the relaxation process of 
the previous paragraph and subtract from the total work 
the energy dissipated as if the system were during the 
whole process in the final stationary state. That is in the 
integral in equation (|15p we replace u(t) with p%. Now 
everything is local. However, even if one could prove that 
a finite quantity is obtained in the limit T — > +00, there 
is no Clausius inequality. In fact, as discussed long ago in 
[ill , a stationary state is not in general a state of minimal 
entropy production. 



In the recent work [17J, Maes and Netocny considered 
the topic of a nonequilibrium renormalized Clausius in- 
equality in the context of a single Brownian particle in 
a time dependent environment. In this setting, by in- 
troducing a potential field such that the corresponding 
stationary state has minimal entropy production, they 
propose a notion of "renormalized heat" which satisfies 
the Clausius inequality. To compare the approach in 17 1 
to the present one, we can consider N independent parti- 
cles in the hydrodynamic scaling limit N — > 00. We then 
obtain a thermodynamic system for which the scheme 
discussed here can be applied. The two notions of renor- 
malization are different. It is not clear to us how the 



approach in 17J can be generalized to cover the case of 
interacting particles in the hydrodynamic scaling limit. 

We are in debt to J. Lebowitz for permanent discus- 
sions on the topics here considered. We also acknowledge 
stimulating exchange of views with C. Di Castro, T. Ko- 
matsu, N. Nakagawa, S. Sasa, and H. Tasaki. 
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